Perturbative renormalization of the first moment of structure functions 

for domain-wall QCD 
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Using the domain- wall formulation of lattice fermions, we have computed the one-loop renormal- 
ization factors of one-link operators which measure the first nontrivial moment of the unpolarized, 
polarized and transversity structure functions, in the flavor nonsinglet sector. The knowledge of 
these factors is necessary in order to extract physical numbers from domain- wall Monte Carlo sim- 
ulations of parton distributions. 

We have automated the perturbative calculations by developing suitable FORM codes. The 
results show that in many instances the total renormalization factors are almost equal to one, and 
that hence the corresponding operators are, for the appropriate values of the Dirac mass M and the 
coupling go, practically unrenormalized. 
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I. INTRODUCTION 



Domain- wall fermions [ESQl provide a solution of the 
Ginsparg- Wilson relation 0] , and as such they possess an 
exact chiral symmetry at nonvanishing lattice spacings 
without at the same time presenting inconvenient fea- 
tures like doublers or nonanaliticities. They constitute 
one of the most promising formulations for simulations 
of chiral fermions on a lattice and for the study of phys- 
ical issues connected with chirality 0. Although Monte 
Carlo simulations of these fermions require more com- 
putational efforts compared with some other nonchiral 
formulations (like Wilson fermions), recently many ad- 
vances have been reported and at present domain-wall 
fermions are widely used in a variety of physical situ- 
ations, for which some of the most recent results and 
investigations can be found in 7, 8., J|, 10^ Jl, J^, Jji] . 

Form factors, structure functions and generalized par- 
ton distributions are also among the phenomenological 
quantities which have been studied by means of simula- 
tions with domain-wall fermions 
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calculation of the perturbative renormalization of the 
operators related to the moments of the deep inelastic 
structure functions, involving the treatment of covariant 
derivatives, has been missing up to now. These renor- 
malization factors, whether perturbatively or nonpertur- 
batively computed, are however necessary for the reli- 
able extractions of physical numbers from Monte Carlo 
simulations of structure functions. The intention of the 
present work is to provide some of these factors from per- 
turbation theory, and we have here considered the lowest 
nontrivial moment of various parton distributions. In 
particular, we present results for the momentum, helic- 
ity and transversity distributions, which give a complete 
description of the quark momentum and spin at leading 



twist. 

This article is organized as follows. In Sect. |n] we 
review the basic perturbative ingredients which are nec- 
essary for the calculations presented in this work, and 
in Sect. IIIII we discuss the peculiar aspects of one- loop 
renormalization for domain-wall fermions. In Sect. IIVI 
then we introduce the operators of which we have evalu- 
ated the renormalization factors, which we give in Sect. 
fVl Finally, in Sect. IVII we make some concluding re- 
marks, and in the Appendix we give the results for the 
quark self-energy and the bilinear operators, where we 
have found some discrepancies with old calculations , 
which however derive only from some constants in di- 
vergent continuum integrals. These discrepancies do not 
affect the expressions of the renormalization factors, for 
which there is complete agreement with Ref. |T^. 



II. PERTURBATIVE DOMAIN- WALL 

We employ the standard formulation of domain- 
wall fermions devised by Shamir 0], where the five- 
dimensional quark action is given by 
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The Wilson parameter is set to r = —1, and the Dirac 
mass M takes values between zero and two (at tree level) 
so that the correct structure of chiral modes (with no 
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doublers) is attained for Ng oo. The chiral projectors 
are P± = (1 ± 75)72. Here and in most of the paper we 
put a = 1. 

The above domain-wall action can be imagined as a 
Wilson action endowed with an additional flavor index 
s plus a special mass matrix for these flavors, explicitly 
given in Eqs. (|10|l to (|13|l further below. The mass ma- 
trix governs the mixing among the flavors and induces a 
sophisticated structure on the flavor space, which at the 
end produces one light quark and TV, — 1 heavy quarks. 
For this reason the tree-level quark propagator turns 
out to have a more complicated form than in the four- 
dimensional Wilson case, and in practical terms pertur- 
bation theory for domain- wall fermions looks like having 
Ng fermion flavors with an involved propagator struc- 
ture in the index s. The gluon flelds and their couplings 
to the quarks are instead kept four-dimensional, that is 
they do not depend on the flfth dimension and are iden- 
tical at each s. The gluon propagator and vertices are 
then just the same as in a four-dimensional lattice the- 
ory. In this work we have used for the pure gauge part 
of the domain-wall action the standard plaquette action, 
and we perform all computations in a general covariant 
gauge, where the gluon propagator is given by 
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where a = 1 and a = correspond to the Feynman 
and Landau gauges respectively. The measure term, the 
gauge-fixing term and the Faddeev-Popov term, as well as 
the quark-gluon interaction vertices, have also the same 
expression as in the Wilson case. Since r = — 1, the ver- 
tices that we need in this paper assume the form (apart 
from color factors) 
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(3) 
(4) 



for the interaction of the quark current with one gluon 
and two gluons respectively, where p stands for the sum 
of the incoming and outgoing quark momenta. 

The construction of the tree-level quar k propagator 
has been started in 0, and [H lH lH H^l and 
then completed and used in in the first full-fiedged calcu- 
lations of the renormalization of the quark self-energy 
and bilinears [isl 123 . .25] . Further perturbative re- 
sults for domain-wall fermions have been obtained in 
in in m El S m nil, and some perturbative cal- 
culations for nonstandard domain-wall actions have been 
carried out in l33l. 



For our perturbative calculations with domain-wall 
fermions we use the same framework of |l8l |25] , where 
one works solely with the dominant contributions to the 
propagators when the number of fiavors Ns (or points 
in the fifth additional dimension) goes to infinity. In this 
limit the decoupling of the chiral modes is exact for mass- 
less quarks, chirality breaking terms being exponentially 
suppressed in the length of the fifth dimension. 

The domain- wall Dirac operator after a Fourier trans- 
form in (four-dimensional) momentum space becomes 

Dst{p) = (5s, t ^47^ sin (5) 

+ {W+ (p) + mM+)P+ + {W-, (p) + mM-,)P^ , 
where the mass matrix is given by 



WfdP) = -W{p)Sg^t + Ss±i,u 



and 



W{p) = 1 -M + 2^sin' 
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In more explicit form. 
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We see that for m = the Af^'s are absent from the 
action and hence the propagator no longer possesses any 
terms directly connecting the two boundaries at s = 1 
and s = Ng. 

In this work we only consider massless quarks. By in- 
verting the above Dirac operator with m = one obtains 
the tree-level quark propagator 
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The expressions of the functions G^{p) and G^{p) are, for large Ns 
A{p) 



(15) 



Gftip) = 7^ f (1 ~ W^(p)e^"(f)) e-(2JV.-.-t)"(p) + (i _ ^(p)e"(p)) e~(^+*)"(p)') + A{p) e-l^-*l"(f) 

F(j)) \ J 
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where a(p) is defined by the positive solution of the equa- 
tion Hill 



and 



cosh(a(p)) 
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2\Wip)\ 



2W{p) sinh(a(p)) ' 
F{p) = l-W{p)e"^P\ 
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These formulae are only valid for positive W, which is 
always the case if < M < 1. When W has a zero, a 
has a logarithmic singularity. For 1 < Af < 2, W can 
become negative if the momentum is small enough. In 
this case the propagator is given by the above equations 
with the replacements 



W 
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which imply that also sinha changes sign. 

To study matrix elements of the chiral modes in per- 
turbation theory, we need to diagonalize the mass ma- 
trix in the fifth dimension. However, since this matrix 
is not hermitian, one has rather to consider the squared 
mass matrix, that is the second-order operators DZ?^ and 
L>tD, which are hermitian and nonnegative and give a 
well-behaved spectrum. In the second-order operators 
the two chiralities are in fact well decoupled. 



The chiral mode is obtained by means of a rotation 
in the fifth dimension of the original quark fields V's(a^) 
to the basis which diagonalizes the mass matrix, and is 
given by 



where 



{x) = ^1 - J2{P+w'o"'M^) + P-W^'-'M^)), 

(22) 



Wo = 14^(0) = 1 - M. 



(23) 



We can see that, because of the damping factors Wq and 
w,^'"", the chiral mode is exponentially localized near the 
boundaries of the fifth dimension. The "physical" quark 
fields that are instead used in Monte Carlo simulations 
are however somewhat simpler, and they are constructed 
only from quark fields exactly located at these bound- 
aries: 



qix) = P+M^) + P^iPnAx) (24) 
qix) = ^j^Ax)P++^,ix)P.. (25) 
We also use these these expressions for our calculations, 
as done in [Ilillll. 

The computation of matrix elements involving states 
and operators constructed from these physical quark 
fields requires additional propagators. We need in fact 
to connect an internal with a physical quark field, and 
the corresponding propagators are given by 



{Ai-P)q{p)) = (e-(^»-^)"(f)p_ +e-(^-i)"(p)p+) ij^sinp,, - e'^^P) (e^^^-i^^^f^P. e-(^=-^)"(f)p+) . (27) 



For our calculations it is also necessary to know their expansions for small momentum 

l-m2 



{<li~P)^^s{p))c = —p- {w^^-'P++wr'P- -J^^^fi{w^'P++w^'^-^P^)), (28) 

{M-pMp))c = {w^^-'P- + w^o-'P+ «-'P- + w^'-'P+) TZ^^^) (29) 
where the factors 1 — Wq are related to the sums of the tree-level exponential damping factors over the fifth dimension: 

JXi / 1 \ 2 1 Finally, we also need the tree-level propagator 

hm V(«;^=-^P++wr^-) =-r^- (30) . ^ 

«=i {Q{-P)q{p)) = — (31) 
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which in the Hmit of small momentum is equal to 



1 — 

{<i{-p)m)c = —r^. (32) 



III. ONE-LOOP RENORMALIZATION 

Matrix elements estimated by importance sampling 
in Monte Carlo lattice simulations need to be properly 
renormalized in order to become meaningful physical 
numbers. They can be considered as (regulated) bare 
quantities, and in order to get physical results they have 
to undergo a lattice renormalization which matches them 
to some continuum scheme. We choose for the continuum 
the MS scheme of dimensional regularization, since com- 
monly Wilson coefficients of operator product expansions 
are computed in this scheme. 

A perturbative lattice renormalization involves both 
lattice and continuum perturbative calculations. At tree 
level, for momenta much lower than the lattice cutoff, 
lattice operators have the same matrix elements as the 
original continuum operators. At one loop one then gets, 
for the case of a multiplicatively renormalized operator, 

{q\0''^'\q) = (l + .g2(-7(0)logaV + i?''^*))-(9|O*''^l9), 

(33) 

{q\d^\q) = (l+5|js(-7(°^log J+i?^))-(9lO"-'=1g), 

^ (34) 
where the lattice and continuum one-loop finite con- 
stants, and i?^^, do not have in general the same 
value, and hence the one-loop renormalization factors 
on the lattice and in the continuum arc in general not 



equal (the one-loop anomalous dimensions are however 
the same). Here and in the following we call for brevity 
f = {gl/lQTT'^)CF (and similarly for g|^), with Cp = 
{Nl - l)/2Nc for the SU{Nc) gauge group. 

The connection between the original lattice numbers 
and the final continuum physical results is given, neglect- 
ing higher-order terms in g^, by M 



(g|0^'|g) 



= 1- 



2 ,,2 



(35) 



where the difference AR = R 
renormalization factor 

Zo{aiJ.,g) = 1 - - 7 



lat 



j^MS determines the 



logaV + Ai? 



(36) 



which converts the lattice operator 0'°* into the physical 
renormalized operator O^^. The computation of these 
renormalization factors requires both lattice and contin- 
uum perturbative techniques (for more details see |35l|'). 
In the domain-wall case it presents additional peculiar 
features that is worth reviewing. 

Let us first consider, in the massless case, the one- 
loop correction to the domain-wall quark propagator 
{q{~p)q{p))c- It can be easily seen that, given the struc- 
ture of the propagators {q{-p)ips{p))c and {ip s{-p)q[p)) c, 
we can write 



{q{-p)q{p))i loop 



wl 



- [i - wD^qipY 



(37) 



(38) 



where 



■^st{p) 



Wo 



)) 



(39) 



■[wq P^+Wq P+- {wq P-+w„^ P+) Y—^ip)- 
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The calculation of the one- loop self-energy diagrams gives 

m 

^stip) = -f [¥{i^P+ + I'P-) + W+P+ + w^P-) ^, 

(40) 

and when the damping factors are also taken into account 
the final result can be written as 

^qip) = 1 ^ 2 (alogaV + Si - ^ Sg) . 

1 — Wq V I — Wq / 

(41) 



Putting all together, we see that the one-loop correc- 
tion to the quark propagator is of the same form as its 
tree-level expression: 



{q{-p)q{p))i loop = 



- {1 - w^)Y,q{p) 



Z-w -^2, 

(42) 



where 

^2 = 1 + (a log a^p^ + Si) (43) 
is the usual quark wave function renormalization factor. 
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whereas 



= 1 - f ^3 = 1+9^ (44) 

is a new feature appearing in domain- wall fermions, 
which represents an additive renormalization to wo, as 
can be seen from 

{l-wl)Z^ ^l^{wo+f^zf+0{f). (45) 

Thus, while the zero mode remains stable under radiative 
corrections, the Dirac mass M = \ — wq \s additively 
renormalized. This effect is due to the terms in 
Eq. 14U|) . which in turn originate from the order a terms 
in the damping factors of Eqs. H28|) and H29|) . We have by 
explicit calculation checked that the part proportional to 
f — a of S3 is zero (its contribution from the half-circle 
diagram exactly canceling the one of the tadpole) , which 
means that E3 and are gauge invariant. The values 
of El, S3 and can be found in the Appendix. 

We remark that in the above domain-wall self-energy 
there is no term Sq proportional to 1/a, which if present 
would signal a breaking of chirality. 

Let us now consider a composite operator q{x) O q(x) 
which is multiplicatively renormalizable. Again, by look- 
ing at the form of the propagators involved, one can see 
that the one-loop matrix element of this operator be- 
tween "physical" quark states is given by 



{{qOq)qq)i loop = 



Aoip)-0 



1 



(46) 



where Ao{p) contains the contribution of the damping 
factors and can be written as 



A, 



oip)^f{- 10 log a'p' + Bo) . (47) 



characterize the quarks in the nucleon: the momentum 
distribution q{x,Q^) (described by the Fi and F2 un- 
polarized structure functions), the helicity distribution 
Aq{x,Q^) (described by the 171 structure function), and 
the (chiral odd) transversity distribution 6q{x,Q^) (de- 
scribed by the hi structure function). They thus provide 
a complete description of quark momentum and spin at 
leading twist. We have also computed the renormaliza- 
tion of the lowest moment of the 32 structure function, 
which receives contributions from twist-3 operators and 
measures the (chiral even) transverse spin. We refer for a 
more detailed discussion of these structure functions and 
in particular of the operators appearing in their operator 
product expansions, some of which are given below, to 
[39. [33 (of which we follow the notation) and references 
therein. 

We have computed the renormalization factors of all 
flavor nonsinglet operators which contain at most one 
covariant derivative. We have chosen in particular 



Ov2,d = q7{iDi}q, 

1 ^ 

Ov2.e = q-fiD^q- - ^qj,D,q, 



(50) 
(51) 



i=l 



which measure the first moment of the momentum dis- 
tributions, 



Oa2,d = qi{il5Di}q, 



(52) 



Oa 



177475 



1 

D4q - o X! ^^^'^sDiq, (53) 



i=l 



which measure the first moment of the helicity distribu- 
tions, 



(54) 



The one-loop expression has the same form as the tree- 
level matrix element. That also the self-energy contri- 
bution to the matrix element fits properly here can be 
seen (for example when = 7^) from inserting Eq. 141|) 
(without S3) in the expression of the contribution of a 
leg in Fig. ^ 



1-wl 1 



1 — Wq 



(alogaV+Si) 



1 - 1 _ ^2 



■7m 



(48) 

which shows that indeed it gives a multiplicative correc- 
tion to the tree-level matrix element: 

1^ . f (a log a V + Si) • 7. • (49) 



which taken together with O^^ determines the first mo- 
ment of the (72 structure function, and finally 



Ot2 = 9Cr4{l75-D2}9, 



(55) 
(56) 



which correspond to the tensor charge and the lowest 
nontrivial moment of the hi transversity structure func- 
tion respectively. We have not explicitly shown the Gell- 
Mann flavor matrices which specialize them to nonsinglet 
operators and hence forbid any mixing with gluonic oper- 
ators, because they are irrelevant for the sake of the cal- 
culation of the renormalization factors. The symbol {} 
denotes symmetrization over the relevant Lorentz indices, 
while [] denotes antisymmetrization. For the covariant 

derivatives D=D — I? we use the lattice discretizations 



IV. STRUCTURE FUNCTION OPERATORS 

The operators that we have considered in this work 
measure the lowest moment of various structure func- 
tions. They include all three parton distributions that 



1 



Df,qix)^- Uf,{x)q{x + fi) - UUx ~ fi)q{x - fi) (57) 



q{x) Df,^ - q{x + fi)Ul{x)~q{x- jl)U^{x- jl) . (58) 
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We have in some cases considered two representatives 
for an operator measuring a given parton distribution. 
They are differentiated and identified by the choice of 
their Lorcntz indices. The lattice operators correspond- 
ing to these choices fall in two different irreducible repre- 
sentations of the hypercubic group (the symmetry group 
of the lattice, the remnant of the Lorentz symmetry), 
and on the lattice they will renormalize in a different 
way (whereas in a continuum scheme their renormaliza- 
tion factors are equal). In Monte Carlo measurements 
one of the two choices can be more convenient to use 
than the other, giving for instance smaller statistical and 
systematic errors, in particular when one considers the 
role played by nonvanishing momenta in numerical sim- 
ulations. 

Since we have done the calculations with Ns — oo, an 
exact chiral symmetry is maintained in all our results, 
and its most important consequence is that the opera- 
tor which measures the lowest moment of the 52 struc- 
ture function does not show any of the power-divergent 
mixings with operators of lower dimension which are in- 
stead present in the case of Wilson fermions. In fact, 
when chiral symmetry is broken Od^ mixes with a lower- 
dimensional operator which in the continuum operator 
product expansion is 

™g97[4757l]g: (59) 

but on the lattice instead has, in place of the mass, a 1/a 
coefficient which becomes infinite in the continuum limit. 
This mixing is forbidden for domain-wall fermions with 
infinite Ns, and O^i is then in this case multiplicatively 
renormalized. In addition, chiral symmetry implies that 
the renormalization constants of corresponding unpolar- 
ized and polarized operators (which differ by a 75 matrix) 
assume the same value. Thus, chiral symmetry gives a 
reduction of the number of independent renormalization 
factors in a given physical situation. 

For operators which contain one covariant derivative 
one needs to perform a Taylor expansion of all vertices 
and propagators at first order in the lattice spacing a 
(which means the external momentum p) . We have cho- 
sen as loop integration momentum the one carried by 
the internal quarks. Choosing it the one carried by the 
gluon would result in much more complicated expressions 
for the order p contributions. 

V. RESULTS 

The diagrams required for the one loop lattice calcu- 
lations of the matrix elements that we have considered 
here are given in Fig. ^ It can be easily seen that all 
tadpole diagrams are diagonal in the fifth dimension, and 
therefore they are equal to the expression calculated with 
Wilson fermions. A leg tadpole has then the value 

Tz=87r2Zo(l-i(l-a)), (60) 




(c) sail (d) operator tadpole 




(c) leg self-energy (f) leg self-energy 

(excluding tadpoles) (excluding tadpoles) 




(g) leg tadpole (li) leg tadpole 



FIG. 1: The diagrams needed for the one-loop renormalization 
of the lattice operators. 

where Zq = 0.154933390231 ... is a well-known integral 
[35|. while the operator tadpoles have the expression 

To = -Ti (61) 

for all operators considered in this work except the tensor 
char ge, for which the operator tadpole vanishes, Tt-^ ~ 
[36L [3^ . The half-circle contribution of the quark self- 
energy as well as the vertex and sail diagrams instead 
have all a nontrivial structure in the fifth dimension, and 
due to their complexity we have deemed necessary to 
compute them using computer programs. 

We have used the algebraic manipulation program 
FORM ^3§\ to construct routines able to carry out all 
needed analytic calculations in an automated way. Sums 
of the damping factors in the fifth dimension and of 
the four-dimensional expressions in momentum space are 
evaluated with Fortran programs. To improve the con- 
vergence of the numerical integrals we use some of the 
techniques given by Liischer and Weisz in |39| . 
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TABLE I: Values for the sums of the vertex and sail diagrams, Vq" ^ , for the momentum, helicity and transversity operators 
considered in this work, in Feynman gauge. The one-loop results of the proper diagrams in a general covariant gauge can be 
inferred from Eqs. ^ to ||7Sj. We remind that V^^;^'^ = V^"j'^ and K^""^ = K^T^- The case of the tensor charge, Vt"^'^, 
has been calculated for the first time in 18]. 
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The one-loop diagrams of Fig. ^ diverge at most log- 
arithmically. Special care is required for the evaluation 
of the divergent terms, for which it is convenient to use 
the method introduced by Kawai et al. in (see also 
[ssf for simple examples). A logarithmically divergent 
integral 

I{p) = j dkl{k,p) (62) 

is expanded in powers of the external momentum and 
split as 

I{p) = J(0) -I- {I{p) - J(0)), (63) 

where 

J(0) = J dkl{k,0) (64) 

is its Taylor expansion to lowest order in p. Since the 
integrals appearing in J do not depend on the external 
momentum, they are much easier to compute on the lat- 
tice than the complete integral of type /. The whole 
dependence on the external momentum remains in / — J, 
which is ultraviolet-finite for a — > and can be evaluated 
by taking the naive continuum limit. Thanks to these 
facts, one is left with computing on the lattice only in- 
tegrals with vanishing momentum, which is technically 
straightforward. It is to be remarked that while / is well 
defined, for finite lattice spacing both J and I— J are sep- 
arately infrared divergent. To compute them one must 
then introduce an intermediate regularization, which we 
choose to be the naive dimensional regularization. The 



associated divergences will at the end cancel out in the 
sum J + (/ — J). 

To summarize, using this method the computation of 
any divergent integral which depends on an external mo- 
mentum is reduced to the computation of lattice integrals 
at zero momentum plus some continuum integrals. 

Identifying and processing divergent terms in an au- 
tomated way for domain-wall fermions turns out to be 
somewhat more complicated and prone to errors than 
for simpler cases like Wilson fermions. We have thus de- 
vised an alternative indirect procedure for the evaluation 
of divergent integrals. This procedure uses the chain of 
equalities 



j-DW _ jDW _j 




_jDW 






J^) 




- J^) 





where DW stands for domain- wall and W for Wil- 
son fermions. In words, we numerically compute the 
difference between the domain-wall and Wilson zero- 
momentum J integrals, and then add the full well-known 
Wilson result (for which several significant digits can be 
obtained without much effort). The key points here are 
that the difference — is a finite lattice inte- 

gral, because the above-mentioned infrared divergences 
exactly cancel, and thus it does not need to be regular- 
ized at all, and moreover that the difference / — J is an 
integral taken in the continuum limit, and so it makes no 
difference whether it is evaluated using domain-wall or 
Wilson fermions. 

We have tested that our procedure works as desired by 
applying it to calculations with overlap fermions, and we 
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have reproduced in this way all results known in the liter- 
ature for the bilinear operators |4lj and the first moment 
of the unpolarized parton distribution [^E 113 ■ This pro- 
cedure is also much more precise than the subtraction 
of a known simple lattice integral with the same diver- 
gent behavior, which has the disadvantage of sometimes 
converging very slowly and hence it requires very large 
integration grids to attain the same accuracy. In a few 
cases we have used the simpler method (which gives just 
a couple of significant digits) for consistency checks. 

The amputated matrix elements that we have calcu- 
lated have on the lattice the form (see Eq. I47|l 



with 



l+5^(-7o log aV + Sc 



Bo = Vo + To + ^i, 



(67) 



(68) 



where Vq is the finite contribution of the vertex and sail 
diagrams (a, b and c in Fig.^, To refers to the tadpole 
arising from the operator (d in Fig. and Ei is the 
finite contribution (proportional to i^) from the quark 
self-energy of one leg, which also includes a leg tadpole 
(e and g, or f and h, in Fig. QJ. We call "proper" con- 
tributions the ones that do not include the self-energy. 
They correspond to the diagrams a-d in Fig. Calling 
1 — a = ^, the one-loop results for them are: 



3 

- e • 6.850272 + T„,,, ) 0^% 



(69) 



2 2 
P 



(70) 



a2:d 



^ 6.850272 -FT,,. d Or.!d 



Opro,er ^ ( | + ^) log a V 



(71) 



V2 ,e 

'5 



e • 6.850272 -fT„,,ejO;,,e: 

Qproper ^ -2 Q j^g ^2^2 (73) 



+K - ^ ■ 6.850272 + Ta„e 



a2 ,e 1 

^proper ^ ^2 ^ ^) ^2^2 (73) 

a = l i- -7 QCn070 I T^, \ ntree 



+Vdr - i ■ 7.850272 + T^, jO^^ 

^proper ^ f (^^^ogaV (74) 

+Vt "=^ - e • 3.792010) c'*['=^ 

(75) 



^proper ^ ^2 + log a V 



a=l 



^ • 6.350272 + rt JO; 



The results due to the sails and vertices, Vq, have been 
for convenience separated in the Feynman gauge values 



Vq"^^, listed in Tabled for various choices of the Dirac 
mass M between M ~ 0.1 and M — 1.9, and the remain- 
ing contributions proportional to ^, which are instead 
independent of M and shown in the above equations. 
Notice that also the tadpoles coming from the operators 
provide contributions proportional to ^. The results for 
the tensor charge, Vt^, are equal to the results for the 
standard tensor current qa^^tl, which were already cal- 
culated in , and are reported, together with the other 
bilinears and the self-energy results, in the Appendix. 

A significant check of our perturbative calculations 
is that the contributions proportional to ^ in covariant 
gauge are constant in M, as already noted in |35ll3^l3^ 
for the case of overlap fermions and discussed more in 
depth in jSJ. Furthermore, they are equal to the re- 
sults obtained with Wilson fermions (and this is the rea- 
son why we can provide more significant digits for these 
contributions). They are also independent of the lat- 
tice representation of the operator (e.g., for Ov2,d and 
Ot,2^e). Their analytic expressions are very complicated 
and highly nonlinear functions of M containing hundreds 
of terms, and the numerical cancellation of this depen- 
dence is a rather strong check on the good behavior of 
the FORM codes, as well as of the integration routines. 

Another reasonably strong check is that the operators 
0^2, d and Oa2,d have the same renormalization constant 
well within the numerical integration errors, as expected 
from chiral symmetry. We have checked that this is also 
true for the pair O^^.e and Oag.e- Thus, in all cases the 
polarized operators have the same renormalization con- 
stants as the corresponding unpolarized operators. Fur- 
thermore, as we have explicitly verified, for the same rea- 
son the 1 /a coefficient of the mixing term of Eq. H59(l aris- 
ing in the one-loop expression of the operator O^i tends 
to zero when the integration grid is refined. This op- 
erator is then for domain-wall fermions multiplicatively 
renormalized, contrary to what happens in the Wilson 
case, where its mixing coefficient goes to infinity in the 
continuum limit. 

In the numerical integration, the convergence can be- 
come slow when M is very close to zero or two. Thus, 
while otherwise a grid of 60 or 80 points in each direc- 
tion is sufficient to obtain about five significant digits, for 
M = 0.1 and M = 1.9 we had sometimes to increase the 
grid to 100 points in each direction in order to achieve 
the same precision. 

To obtain the complete one-loop amplitudes we have 
now to add to the results of the proper diagrams the 
1-loop contributions of the self-energy which are propor- 
tional to iy. 



Si = 5' ( (1 - log a V + + f • 4.792010 ) , (76) 



where 



10.8750 when M = I, while for other 
values of M the Feynman-gauge finite terms Ej^"^^ are 
given in the Appendix. The complete one-loop lattice 
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TABLE II: Values of the renormalization factors for the various operators, for domain- wall QCD at = 1 and /j, = 1/a. We 
remind that Za2,d — Zv2,d and Za2,e = Zv2,e- 



M 


Zv2,d 


Zv2 


Zdi 


Zt, 


Zt2 


0.1 


0.9979 


0.9945 


0.1931 


0.8494 


9952 


0.2 


0.9984 


0.9948 


0.6480 


8529 


9963 


0.3 


9988 


0.9951 


0.7958 


0.8557 


0.9972 


0.4 


n 9992 


n 9952 


8683 


8582 


9980 


0.5 


n 999 


995S 


0.9112 


8604. 


9986 


0.6 


0.9997 


n 9953 


9396 


0.8624 


0.9991 


0.7 


9998 


0.9951 


9600 


0.8643 


9995 


0.8 


9999 


0.9949 


0.9755 


0.8661 


9998 


0.9 


n 9998 


0.9946 


0.9879 


0.8678 


1.0000 


1.0 


0.9997 


0.9941 


0.9984 


0.8694 


1.0000 


1.1 


0.9994 


0.9935 


1.0077 


0.8710 


1.0000 


1.2 


0.9990 


0.9927 


1.0164 


0.8725 


0.9998 


1.3 


0.9984 


0.9917 


1.0250 


0.8740 


0.9995 


1.4 


0.9976 


0.9905 


1.0342 


0.8754 


0.9990 


1.5 


0.9966 


0.9890 


1.0446 


0.8768 


0.9982 


1.6 


0.9952 


0.9872 


1.0578 


0.8782 


0.9973 


1.7 


0.9934 


0.9849 


1.0765 


0.8797 


0.9959 


1.8 


0.9911 


0.9819 


1.1086 


0.8813 


0.9942 


1.9 


0.9878 


0.9780 


1.1900 


0.8832 


0.9917 



results are then, for M — 1: 



'^V2.d 



o: 



,lat 
a2 ,d 



o 



lat 



o 



a2 .e 
lat 



Qlat ^ 



^+f{l log a V - 4.4059 + i) ) 0%% (77) 
l+f{^-\og a V - 4.4059 + i) ) (78) 
l + log aV - 3.7445 + ) 01^1. (79) 
l+f{^ log aV - 3.7445 + ) Of^te , (80) 
l + g2•0.1845)o^7^ (81) 



1 + 5' ( log a V + 14.4633 + ) O^' - (82) 
l + 5'(31ogaV 



5.0052+ -^))0*n83) 



To establish the connection with the corresponding con- 
tinuum quantities we also need to know the one-loop am- 
plitudes for the same operators in the MS scheme [SalST^ : 



^V2 



= 1+.9 



, /», 40 ^ 



r\tree 
^V2 ' 



MS 



o. 



Otree 
di ' 



MS 



^t2 



l+^^(31og^-5+^,))o*: 



(84) 

(85) 
(86) 
(87) 

(88) 



Putting all together, we obtain the factors that allow the 
matching from the domain-wall lattice theory to the MS 



continuum scheme, for M = 1: 

l-.92(^logaV 
1 



^V2.d 



nMS 

^a2,d 



O 



MS 



o 



MS 



MS 



o, 



QMS ^ 



nMS 

^t2 



^(flogaV^ 
^(^logaV^ 
-(^logaV^ 



0.0386jjO^,"^^„ (89) 
0.0386)) Oi-^^^,, (90) 
0.6999)) 0^f„ (91) 



0.6999jjOi",%, (92) 
, ... (93) 
1 - .g^^logaV^ + 15.4633)) Ojf, (94) 



f ■ 0.1845 lOf,"*, 



1 - .g^(31ogaV^ - 0.0052) )0: 



-\lat 



(95) 



Notice that the part proportional to ^ has canceled be- 
tween the lattice and continuum expressions, and these 
renormalization factors are hence gauge invariant. 

For simulations of domain- wall QCD at go = Ij setting 
/i = 1/a one obtains the values 



r,MS 

^V2M 



^MS 



= 0.9997-0: 



lat {domain —wall, Af— 1.0) 
V2 ,d 



(96) 



= 0.98920 -O: 



lat (Wilson) 
V2 ,d ' 



(97) 



O 



MS 



o 



MS 



Q 0007 (domain— wall, Af— 1.0) 

a2~d 

^ 0.99709 • 0"^*i'^^'^°"\ 

a2 ,a ' 
_ 9941 0^^^ (domain— wall, M— 1.0) 

= 0.97837 • 0^'^^y^''"°"\ 

Q QQ^^ 0lat (domain— wall, M— 1.0) ^99) 

= 0.99859 • Ol''/ J^*'"°"\ 
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qMS _ Q ggg^ Q^"-* (domain-wall, M=1.0) (IQO) 
qMS _ Q ggg4 . ^'at (domain-wall, M=1.0) (IQl) 

= 0.85631 • O^";* ('^^'-"), 

qMS _ ^ 0000 • O'"* (domain-wall, M=1.0) (102) 

= 0.99559 -Ojf (^^'^°"), 

where for comparison the corresponding Wilson results 
are also shown. Of course the domain-wall renormaliza- 
tion factors vary with M. For example, for M — 1.8 
(which is almost at the edge of the allowed parameter 
space) their values are instead 





= 0.9911 


^lat {doinain 


—wall. 


M= 


1.8) 


(103) 


nMS 

'^a2,d 


= 0.9911 


^lat {domain 
^a2,d 


-wall, M= 


1.8) 


(104) 


qMS 

V2 ,e 


= 0.9819 


r^lat {domain 


— wall. 


M= 


1.8) 


(105) 


a2 ,e 


= 0.9819 


r^lat {domain 
a2,e 


—wall. 


Af= 


1.8) 


(106) 


nMS 


= 1.1086 


^lat {domain 
^di 


—wall. 


M= 


1.8) 


(107) 


/-)MS 


= 0.8813 


y^lat {domain 


—wall. 


Af= 


1.8) 


(108) 


^MS 

"^t2 


0.9942 


^lat (domain 

^t2 


—wall. 


Af= 


1.8) 


(109) 



Results for other choices of M and 170 can be easily ob- 
tained from the numbers presented in this Section. Table 
nil contains the values of the renormalization factors as a 
function of M for go = 1- 

All the renormalization constants presented in this 
work can also be used in unquenched simulations, pro- 
vided that one computes only matrix elements of flavor 
nonsinglet quark operators, for which at one-loop internal 
quark loops never appear. The numbers for the transver- 
sity operators can be however considered unquenched 
even in the singlet case, since no chiral-odd gluon op- 
erators exist, which would constitute the only possibility 
for having a mixing. 

We can easily notice that the renormalization correc- 
tions that we have obtained for domain- wall fermions are 
in general small. In particular when the Dirac mass is 
M = 1 or not too far from it they are not too different 
from the corresponding Wilson results. In many cases the 
total renormalization factors turn out to be quite close 
to one, reflecting the fact that the domain-wall one-loop 
amplitudes have almost the same values as the corre- 
sponding MS results. The only exception is the tensor 
charge, which on the other hand is also the only case 
which for fermions which break chiral symmetry cannot 
be renormalized, because of its power-divergent mixing. 
From this point of view, domain- wall fermions appear to 
behave quite at variance with overlap fermions, for which 
the renormalization factors are generally not small, giv- 
ing in many cases rather large one-loop corrections to 
the tree-level matrix elements [s^ Is^, 0|. The origin 
of most of these large effects can be traced back to the 
El contribution from the half-circle diagram of the self- 
energy, which for overlap fermions is rather big. On the 



contrary, for domain-wall fermions Si does not deviate 
too much from the Si of Wilson fermions. Thus, apart 
from Odi when M is away from one, the renormalization 
factors computed in this work give small corrections at 
the standard accessible couplings. 

We have computed the bilinears and the self-energy 
anew, and we have found some discrepancies when com- 
paring our results, which we report in the Appendix, with 
the numbers given in [lH. These discrepancies derive 
only from the continuum integrals that are needed to 
compute the divergent terms with the Kawai method (for 
S3 and the tensor charge, which are finite, we completely 
agree). The differences between our results for the i?'"* 
quantities and those in |l8l | are indeed in all cases inde- 
pendent of M and always are an integer or half-integer 
number. The MS renormalization factors given in [13 
are also different, and for the same amount, from what 
is found elsewhere in the hterature (e.g., [s^ |33, 0|)- 
All these differences cancel then in the expressions of the 
renormalization factors, which do not present any dis- 
crepancies with ^3 ■ 

VI. CONCLUSIONS 

In this paper we have presented the computation of the 
one-loop renormalization factors of a few operators which 
measure the first nontrivial moment of various structure 
functions, giving a complete description of the quark mo- 
mentum and spin at leading twist. We have used domain- 
wall fermions, and the associated chiral symmetry plays 
an important role in the structure of the strong radiative 
corrections. 

We have automated the calculations by developing 
suitable FORM codes. This will make it easier to con- 
sider the case of more complicated operators. The renor- 
malization factors that we have found turn out to be in 
many cases close to one. 
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APPENDIX: SELF-ENERGY AND BILINEARS 

We report here the results for the quark self-energy 
and the bilinear operators, which were first calculated in 
Ts^ in Feynman gauge. 

Taking into account that the results for the pseu- 
doscalar and axial- vector operators are equal to the ones 
for the scalar and vector operators respectively, and that 
the tensor operator has been reported in the main body 
of the paper, the one-loop results for the proper diagrams 
that we need here are: 
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TABLE III: Values of the domain-wall constants needed for the renormalization of the self-energy and the bilinear operators. 



We remind that Vp"=i = Vg"^^ 


and = Vv > s-nd that 


0^ V //I ^ 2\ 






M 




S3 


Zw 






0.1 


11.6603 


51.0482 


-483.6145 


7.8219 


6.3355 


0.2 


11.5099 


50.7450 


-225.5333 


8.6070 


6.3380 


0.3 


11.3829 


50.4885 


-100 K K 

-138.5959 


r\ C\ A C\ A 

9.2424 


6.3408 


0.4 


11.2730 


50.2664 


-94.2495 


9.8020 


6.3438 


0.5 


11.1772 


50.0726 


-66.7635 


10.3176 


6.3472 


0.6 


11.0939 


49.9038 


-47.5274 


10.8074 


6.3509 


0.7 


11.0221 


49.7582 


-32.8076 


11.2834 


6.3549 


0.8 


10.9616 


49.6352 


-20.6813 


11.7549 


6.3594 


0.9 


10.9124 


49.5351 


-10.0071 


12.2296 


6.3644 


1.0 


10.8750 


49.4588 


0.0000 


12.7151 


6.3699 


1.1 


10.8504 


A C\ A r\C> A 

49.4084 


9.9815 


13.2189 


6.3762 


1.2 


10.8399 


49.3865 


20.5777 


13.7496 


6.3831 


1.3 


10.8455 


49.3972 


32.5696 


14.3176 


6.3910 


1.4 


10.8699 


49.4461 


47.0915 


14.9360 


6.4000 


1.5 


10.9170 


49.5411 


66.0549 


15.6225 


6.4102 


1.6 


10.9923 


49.6935 


93.1753 


16.4024 


6.4219 


1.7 


11.1041 


49.9198 


137.0349 


17.3150 


6.4356 


1.8 


11.2652 


50.2462 


223.3164 


18.4274 


6.4516 


1.9 


11.4979 


50.7176 


480.4824 


19.8800 


6.4706 



0proper -2 



2^2 



(-4 + logaV 



(A.l) and thus 



-hFs"=^ - C ■ 5.792010 )0^'^*=^ 



^proper -2 



.92((-l + 01ogaV (A.2) 
- ^-4.792010) 0( 



^tree 

W ■ 



There is no operator tadpole for the bilinears, and adding 
the self-energy contribution proportional to (Eq. (|76|) ') 

El =52 (^(1-0 log aV + Si"=i+C- 4.792010), (A.3) 
we get, for M = 1, 

O^f = (}-f{^ log a V - 23.5901 + ^) ) O^lh-^) 
O^f = {l^f- 17.2450)of;*=^ (A.5) 

The one- loop results in the MS scheme are [s^ [s^, : 

OT = (l+.9'(-31og^ + 5-e))or^(A.6) 



OMS r^tree 
V — ■> 



(A.7) 



0.8430 • O 



lat {domain— wall, 1.0) 



0.89064- O' 



lat {Wilson) 



o 



MS 



0.8430 • O 



lat {do'main — wall, 1.0) 



(A., 



(A.9) 



0.80922- O 



lat {Wilson) 



qMS _ Q . Q^^^ {domain-wall, M^l.O) ^q-j 

= 0^82592-0;;* (w-^i.on)^ 



Oo/^/^/^o y^lat (Wilson) 
.86663-6'^ ^ 

where for comparison the corresponding Wilson results 
are also shown. The renormalization factors for other 
values of M and go can be easily obtained from the num- 
bers given in Table Uni where we in addition to Si we 
report also the results for the quantities S3 and 



2wo 



(A.12) 



which determines the amount of additive renormalization 
to Wo = 1- M. 
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